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Abstract 

We have carried out a NLO analysis of the world data on polarized DIS in the 



MS scheme. We have studied two models of the parametrizations of the input 
parton densities, the first due to Brodsky, Burkhardt and Schmidt (BBS) which 
gives a simultaneous parametrization for both the polarized and unpolarized 
densities and in which the counting rules are strictly imposed, the second in which 
the input polarized densities are written in terms of the unpolarized ones in the 
generic form /S.q{x) = f{x)q{x) with f{x) some simple smooth function. In both 
cases a good fit to the polarized data is achieved. As expected the polarized data 
do not allow a precise determination of the polarized gluon density. Concerning 
the polarized sea-quark densities, these are fairly well determined in the BBS 
model because of the interplay of polarized and unpolarized data, whereas in the 
second model, where only the polarized data are relevant, the polarized sea-quark 
densities are largely undetermined. 
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1. Introduction. 

Experiments on polarized deep inelastic scattering (DIS) were initiated at SLAC 
by the SLAC- Yale group |]T| soon after the discovery of Bjorken scaling. Enormous 
impetus was given to the subject by the European Muon Collaboration (EMC) exper- 
iment at CERN in 1988 whose results seemed to imply a "spin crisis in the parton 
model" 0. Much theoretical and experimental work has followed and today there is a 
rich program of experiments under way (E154, E155 at SLAC; HERMES at HERA)or 
in the progress of being set up (COMPASS at CERN). 

Experiments on unpolarized DIS provide information on the unpolarized quark den- 
sities q{x, Q^) and gluon density G{x, Q"^) inside a nucleon. Polarized DIS experiments, 
using a longitudinally polarized target, give us more detailed information, namely the 
number densities of quarks q{x, Q^)± and gluons G{x, Q'^)± whose helicity is re- 
spectively along or opposite to the helicity of the parent nucleon. The usual densities 
are 

q{x, Q2) = g^(x, + g2) ^ ^^^^ ^ ^^(^^ g2) ^ ^ (^^ g2) 

and the new information is then contained in the polarized structure function gi{x, Q^) 
which is expressed in terms of the polarized parton densities 

Aq{x,Q^) = q+{x,Q^)-q4x,Q^) , AG{x,Q^) = G+{x,Q'') - G4x,Q^) . (2) 

Two developments in the past few years have made it possible and worthwhile to 
attempt a detailed comparative study of the polarized parton distributions. On the 
one hand, a wealth of new data, much of it high quality, has appeared [4 - 12]. On the 
other, the theoretical calculation of the Altarelli-Parisi splitting functions to two-loop 
order has been, after several hiccoughs, completed successfully [pT^ . 

It would be wrong, however, to imagine that the polarized densities can now be 
determined to the same accuracy with which the unpolarized densities are known. 
This can be understood quite simply. Up to the present the polarized data consist 
solely of fully inclusive neutral current (in effect, photon induced) reactions on protons 
and (via deuterium or Helium-3) on neutrons, i.e. one has information on the two 
polarized structure functions gi{x, Q"^) and 5'"(x, Q"^) . Even if one makes some 
simplifying assumptions about the polarized sea, e.g. A'u(x) = /S.d{x) = As{x) 
or Au{x) = Ad{x) = 2As{x) one is still expressing two experimental functions in 
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terms of four densities: Au{x,Q'^), Ad{x,Q'^), Aq{x,Q'^) and AG{x,Q'^). What is 
lacking here is the information from charged current reactions which plays an important 
role in pinning down the unpolarized densities. The situation is alleviated by the 
beautiful connection between the first moments of the polarized parton densities and 
weak interaction physics. Namely, one has the connection with neutron f3 -decay, via 
the Bjorken sum rule, 

dx\Au{x, Q^) + Au{x, - Ad{x, Q^) - Ad{x, Q^)] = gA/gv (3) 

JO 

and, to the extent that flavour SU{?)) is a good symmetry, the connection with 
hyperon (3 -decay 

C dx{Au{x,Q^)+ Au{x,Q^) + Ad{x,Q^)+ Ad{x,Q^) 
Jo 

- 2[As{x,Q^) + As{x,Q^)]} = 3F-D. (4) 



The values of gA/gv and 3F-D are taken from [jT4| 



gj^/gy = 1.2573 ± 0.0028, 3F - D = 0.579 ± 0.025 . (5) 

Assuming a roughly flavour-independent polarized sea, allows one to interpret Eqs. 
d^) and as statements about the first moments of the polarized valence densities: 

dx[Au^{x, Q2) - Ad^{x, Q^)] = 1.26 , (6) 

JO 

r dx\Au^{x, Q^) + Arf^(x, g^)] = 0.58 , (7) 

JO 

which immediately suggests that A(i„(x,Q^) is of opposite sign to Au^^x^Q"^) and 
of roughly comparable magnitude, in agreement with simple SU{Q) models of the 
proton wave-function. Eqs. and (|) are crucial supplements to the polarized DIS 
data. 



In seeking input parametrizations of the polarized densities into the QCD evolution 
equations one must clearly respect the positivity of the number densities 5'(a;,(5o)± , 
which via (|l]) and (0) is equivalent to demanding 

\Aq{x,Ql)\<q{x,Ql) (8) 

There are, in addition, certain counting rules relating to the behaviour of Aq{x) / q{x) 
as a; ^ and x —>■ 1 , which follow in the parton model from perturbative QCD and 
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the form of the infinite momentum frame nucleon wave function [15 - 17]. 

We have examined two classes of models for the input densities. The first, due 
to Brodsky, Burkardt and Schmidt [0 is unusual since it directly parameterizes 
Q^) and g_(x, Q^) at some = (rather than Ag(x, Qo) and q{x,Ql)) 
so that the free parameters are determined from a simultaneous fit to the polarized 
and unpolarized data. Positivity is simple to implement and the counting rules are 
imposed exactly at Q'^ = Ql ■ 

In the second, each polarized parton distribution is written in the generic form 
Aq{x,Ql) = f{x)q{x,Ql) at some Ql , with the usual densities q{x,Ql) determined 
from the unpolarized data (in practice we utilize the MRS(A ') set of distributions p!8[]). 
The functions /(x) are parameterized so as to respect positivity but the counting rules 
are not imposed in a strict fashion. 

Finally it should be remembered that beyond the leading order in perturbation 
theory the parton densities become scheme dependent. In this paper we work in the 
MS scheme. 

In Section 2 we wish to draw the reader's attention to certain interesting qualitative 
features of the polarized DIS data. In Section 3 we explain the method of analysis and 
in Section 4 discuss the parametrization of the models and their properties. Our results 
are presented and analyzed in Section 5 and conclusions follow in Section 6. 

2. Implications of qualitative features of the data 

The structure functions gi'"'{x, Q"^) are determined in the following x range: 0.003 < 
X < 0.8 . For small x there are large errors on the measured values, but if one takes the 
central values of the data points as indication of the trend of the behaviour as x ^ 
then one is led to some surprising conclusions. To see this, consider the expressions for 
the structure functions gi'"'{x, Q"^) . They can be expressed in terms of contributions 
Ag3(x,Q^), Ag8(x, Q^) and AS(x, Q^) of definite fiavour symmetry as 

gr{x, Q') = {±^Aq,{x, Q') + ^Aqsix, Q') + ^AJ^ix, Q')}{1 + 0(a,)} (9) 
12 6b 9 

where the fiavour non-singlet contributions are 

Aqsix, Q^) = Au{x, g^) + Au{x, Q^) - Ad{x, Q^) - Ad{x, Q^) , (10) 

Aqs{x,Q^) = Au{x,Q^) + Au{x,Q^) + Ad{x,Q^) + Ad{x,Q^) 

- 2[As{x,Q^)+As{x,Q^)] (11) 
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and the singlet contribution is 



AS(x, Q') = Q') + Aqfix, Q')] . (12) 

/ 

The gluon contribution to gi{x,Q'^) is hidden in the 0{as) correction terms in 
(H) and a more precise expression will be given in Eq. (|T8l). 

One sees then that the difference gi{x,Q'^) — g^i^x^Q"^) is a purely non-singlet, 
whereas the sum gi{x, Q"^) + (7" (x, Q^) is a mixture of singlet and non-singlet contri- 
butions. 



Now according either to the small-x behaviour of the evolution equations 19 



or 



to the summation of double logarithmic terms at small x |2^, the singlet contribution 
should dominate over the non-singlet terms as x — , which, by the above, could 
imply ( gi + g"^) > {g{ — g"^) at small x. The data (see Fig. 1) show precisely 
the opposite trend. However, while the theoretical arguments predict the form of the 
behaviour as x —>■ , namely CnsX'"'"" or CsX~""' with > a„s > 0, the values of 
the coefficients Cns, Cs are sensitive to the structure of the parton distributions at 
Qq. So there need not be a contradiction at presently measured x- values. However, as 
experiments probe smaller and smaller x, there ought to be a dramatic change in the 
trend of the data or else the theoretical arguments are incomplete. 

On a practical level concerning the present data, the behaviour in Fig. 1 can lead 
to difficulties when one makes, as one is forced to do, some simplifying assumption 
for the polarized sea, such as Au{x,Ql) = Ad{x,Ql). For then Aq3{x,Q^) in ([lOD 
practically depends only on valence distributions and in some cases one finds that the 
best-fit parameters tend to make Ady{x,Q'^) so large and negative at small x that 
positivity is violated. 

3. Method of analysis 

Measurements of polarized deep inelastic lepton nucleon scattering yield direct in- 
formation on the virtual photon- nucleon asymmetry A^{x, Q^) . Neglecting as usual 
the subdominant contributions, A^{x,Q'^) can be expressed via the polarized struc- 
ture function g^{x,Q^) as 



Q') = ^wrfhk = ^wB^lMi + R'^ix, Q')] , (13) 

where 



j^iy ^ _ 2xF^)/2xFi^ (14) 
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and and are the unpolarized structure functions. 

Usually the theoretical analysis of the data is presented in terms of gi{x^Q^) 
extracted from the measured values of Q^) according to ([T3|) using different 

parametrizations of the experimental data for F2 and R. The parametrizations |^ of 



F2 and p2| of R were used in the most recent analyses. Some experimental groups 
assume also scaling of A^{x,Q'^) in their extraction of g^{x,Q'^) . However, 
bearing in mind the recent NLO calculations of gi in QCD, this assumption is not 
theoretically correct, especially in the small x region. 



In our analysis we follow the approach first used in |2^, in which the next-to- leading 
(NLO) QCD predictions for the spin-asymmetry A^{x,Q'^) are confronted with the 
directly measured values of {x, Q"^) rather than with the gi{x, Q"^) derived by the 
procedure mentioned above. A further advantage of such an approach is that higher 
twist contributions are expected to partly cancel in the ratio ([13D , in contrast to the 
situation for gi{x, Q^) . Usually to avoid the influence of higher twist effects Q^- 
cuts in the (x, Q^) data set are introduced. Bearing in mind that in polarized DIS 
most of the small x data points are at low , a lower than usual cut is needed 
{ > 1 GeV^ ) in order to have enough data for a theoretical analysis. We consider 
that in this approach such a low Q^-cut is more reasonable. 

In NLO approximation 

4iV/ ^ 2xg^{x,Q^)NLO ..^x 

In (|15D N = p, n and d= {p + n)/2 . 

To calculate gi{x,Q'^) and 2xF^{x,Q'^) we have used the analytic NLO solutions 
for the moments in Mellin space with the nth moment being defined by 

M^{Q^) = r dxx^'-^xg^ix^Q^) , n = l,2,... (16) 
Jo 

(Q2) = dxx'^-^2xFi{x, g2) , n = 2, 3, ... (17) 

JO 

In NLO approximation with = 3 active flavours the moments ([16|) of the 
structure function gi{x, Q^) can be written in the form [^]: 



+ ^5C„^AG(n,g2)-,o" , ( 
/vr 
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MtiQ') = \[M^r^\Q') + Ml"Hg')](l - . (19) 

where Aga, Agg and AS are the moments of the flavour non-singlet and singlet 
combinations (|10|), ( ]TT| ) and (]T2|), respectively, while Q^) denotes the moments 

of the gluon density (j(a;,(5^). The subscript "LO" in (0) means that the moments 
of the corresponding densities satisfy the LO Q^-evolution equations, in which for the 
strong coupling constant as{Q'^) the NLO approximation 

Air f3olnQyA^ (3l InQ^A^ ^ ^ 

is taken. In ([T9|) for the probability of the deuteron to be in D-state we have taken 
LJd = 0.05 ± 0.01, which covers most of the published values . 

All quantities - the anomalous dimensions Sj^j up to two-loop approximation and 
the moments of the coefficient functions 6C!^ and 6C^ in one-loop approximation 
needed to derive the analytic "LO" and NLO solutions for the moments of the parton 
densities, can be found, for instance, in |p3i. 



Unlike paper |^ where the expressions for the moments of gi and 2xFi are 
numerically Mellin-inverted to yield the structure functions in Bjorken x-space, we 



follow a method which presents the structure functions analytically. Having 

the NLO Q^-evolution of the moments (|18]) and (|19[) we can write the structure function 
gi{x,Q^)NLO in the form: 



Nrr 



j=0 



n=0 



(21) 



where Q"^{x) is a set of Jacobi polynomials and Cj (a,P) are coefficients of the 
series of G"'^(x) in powers in x: 

n 

Q"J{x)=Y.cf\a,P)x^. (22) 

Nmax, « and /? have to be chosen so as to achieve the fastest convergence of the 
series on the R.H.S. of Eq. (pi]) and to reconstruct xgi with the required accuracy. 
We use a = 0.7, /3 = 3.0 and Nmax = 8 . These values guarantee an accuracy 
better than 2.5% in the experimental x range: 0.01 < x < 0.8 and better than 5% 
for smaller x. 
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The same method has been apphed to calculate the unpolarized structure functions 
2xFf{x,Q^)NLO from their moments (|l^. Following the results of Ref. we use 
for the quantities N^ax, « and /? in this case: 

N^ax = 12, a = -0.85, f3 = 3.0 

in order to guarantee an accuracy better than 10~^ in the x range mentioned above. 
In the present calculations of 2xFi the MRS (A') parametrization for the input un- 
polarized parton densities has been used. 



As already mentioned in the introduction all calculations are performed in MS 
scheme. In this scheme the first moment of the Wilson coefficient function 6Cf = , so 
the axial anomaly contribution of AG'(1, Q"^) to the first moment [Q"^) = T^IQ"^) 
of gi{x,Q'^) vanishes, z.e. the singlet axial charge ao((5^) coincides with AS(1,(5^) 
(twice the total helicity) carried by the quarks in the nucleon. In the present study we 
use the MS scheme because of our choice of the input parton distributions (see, Eqs. 
( [3T| ) and (0)). A NLO analysis of the polarized DIS data in a renormalization scheme. 



which allows for an anomalous gluonic contribution AG(1, Q ) to the Ellis- Jaffe sum 



rule [53] is presented in p 



The value of A;^ is taken to be Arj^ijif = 3) = 284 MeV. This value corresponds 



to ^jisin-f = 4) = 231 MeV |]I8[ according to the requirement be continuous 



across each threshold > m^. The contributions of charmed quarks to gi{x,Q'^ 



is assumed to be neghgible at present energies ||30| and will not be considered in our 
analysis. 

The last step before fitting the theoretical predictions for Ai{x, Q"^) to the data is 

2 

' 



to choose the input polarized parton densities at some fixed value of = Ql , evolve 



them to Q and then put them into (pTI). 

4. Models for the input parton distributions 

We have studies two classes of models: 

(a) the model of Brodsky, Burkardt and Schmidt (BBS) |T3, which directly param- 
eterizes the parton densities q±{x) and which respects the perturbative QCD counting 
rules exactly 

(b) an example of parametrization of the form Aq{x) = f{x)q{x) where q{x) are 



the unpolarized parton densities of Martin, Roberts and Stirling MRS(A') |jT8[ and in 
which the counting rules are somewhat relaxed. 
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(a) The BBS counting rule model 



It is well known that the valence Fock states with the minimum numbers of con- 
stituents dominate in determining the behaviour of the unpolarized valence distribu- 
tions dy{x) in the region x — > 1 and this leads l3l[], via perturbative QCD, to 
the prediction 

Uy{x), d^{x) ^ {I- xf (a; ^ 1). (23) 



The same arguments applied to the polarized case suggest |jT5|, jl^ that 

g,^(x) ~ (l-x)^ g„_(x) ~ (1-x)^ (x ^ 1) (24) 
implying, via (^31) , that 

(x-l). (25) 

Simple perturbative arguments based upon the splitting functions for q — > qG 
and G ^ qq can then be used to predict the behaviour of the gluon and sea densities 
generated from the valence quarks in the region x ^ 1 , namely |jT6|, |17 



G(x)~(l-x)^ ^1^-1 (^-1) (26) 

and 

qix) ~ (1 - x)\ ^ ^ 1 ix-^l) . (27) 

q 

Unfortunately these simple sum rules are not compatible with the evolution equa- 
tions. If they hold at some Q"^ then the power of (1 — x) involved will grow like 
InlnQ^ at large . It is therefore not clear to what extent the above rules should 
hold at some at which one chooses to parameterize the parton densities. Presum- 
ably they should be more accurate at momentum scales Qq ~ ^qcdi so that for the 
typical values of Ql utilized in the data analyses, 1 < — ^ GeV"^, one might expect 
somewhat higher powers of (1 — x) to appear. Indeed, analyses of the unpolarized 
data yield typically, at Qq = 4 GeV^, 



(l-x)(3-^), 4(x)~(l-x)(^-^) (28) 



somewhat at variance with (^) and with the faster decrease of dv{x) being related 
to the experimentally established behaviour of -F^(x)/-F2 (x) as x ^ 1. 

The BBS model is interesting because it asks how well the data can be fitted if the 
exact conditions (|23| - ^) are imposed. 
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Concerning the behaviour as a; — > 0, perturbative arguments based on the sphtting 
functions predict fl^ that the gluons and antiquarks 

Aq(x) AG(x) 

^tV ' ^vV ^ ^^^^^-^ (x 0) . (29) 
q[x) G[x) 

This behaviour is consistent with Regge-type arguments. 

For the gluons BBS use wave-function arguments to suggest that the const in ( |29D is 
approximately equal to 1 and their gluon density automatically satisfies this constraint. 
For the sea the const is determined by the fit to the data. 



In the original BBS paper [|l^] an analysis of the polarized data then available and a 
fit to the MRS (DO ') parametrization of the unpolarized data were performed, but with- 
out taking account of the corrections from QCD evolution. As on the one hand, there 
is now much more polarized data available, especially at smaller x, and the MRS(DO') 
fit has been superseded by other parametrizations which fit the unpolarized low-x data 
far better, and on the other, the NLO calculations in QCD have been completed, we 
have felt it important to examine to what extent the original claims of the BBS model 
remain valid. 

The input helicity-dependent parton densities at Qo in the BBS model have the 
following form: 

x{Au{x) + Au{x)) = x^-''^{l-xf[A^ + B^{l-x)-Cu{l-xf -D^{l-xf] , 
x{Ad{x) + Ad{x)) = x^-'"'{l~xf[Ad + Bd{l-x)-Cd{l-xf-Dd{l-x)^], 
x{As{x) + As{x)) = x^~'''{l-xf[As + Bs{l-x)-Cs{l-xy -Ds{l~xf] , 

xAG{x) = x^-'''{l-x)^[Ag + Bg{l-x)][l-{l-x)^] , (30) 

while the unpolarized parton densities at Qq are given as 

x{u{x)+u{x)) = x^-'^^{l-xf[Au + Buil-x) + Cuil-xY + Duil-xf^_ 

x{d{x) + d{x)) = x'-"''{l-xf[Ad + Bd{l~x)+Cd{l-xy + Dd{l-xf] 

x{s{x) + s{x)) = x^-"''{l~x)^[As + Bs{l-x) + Cs{l-x)^ + Ds{l-xf^ 

xG{x) = x'~''^{l-x)''[Ag + Bg{l-x)][l + {l-xY] . (31) 

The constraints 

Ag + B, = Gg + D, (32) 

are imposed on the constants Ag, Bg, Gg and Dg in (|30D and ( plj ) in order to ensure 
the convergence of the helicity-dependent sum rules (^ and (^). Thus in the BBS model 



5 
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the Regge behaviour of the polarized quark densities Ag ~ a;~°^ is automatically one 
unit less than the unpolarized intercept aq : aR = aq — 1. Isospin symmetry at low x 
requires 

A^ + B^ + Cu + D^ = Ad + Bd + Cd + Dd. (33) 

If in addition to (|32|) and (|33|) the helicity- dependent sum rules @ and (§) and the 
energy-momentum sum rule for the unpolarized densities (|3TD are taken into account, 
the number of the unknown parameters associated with the input polarized densities 
is reduced to = 16 — 7 = 9. These free parameters - we have chosen 

{Bu, Cu, Du, Cd, Cs, Ds, Bg, ttq, ttg} , (34) 

are determined by a simultaneous fit of the theoretical predictions (^) and the BBS 
unpolarized parton densities ( ^T]) to the world A^(x, Q^) data and the MRS(A') set 
of unpolarized parton densities, respectively. We recall that the MRS (A') densities at 
next-to-leading order are parameterized at Qq = 4 GeV^ and are well determined 
from the wide set of unpolarized experiments. Note that the choice of unpolarized 
densities is not crucial for the present analysis. 

(b) Models based on the unpolarized parton densities 

We have also analyzed the polarized DIS data using the following expressions for 
the input polarized parton distributions at 

x/S.u^(^x, Qq) = Tjy^Ay^x ^xu^{x^Qq) , 

xM^{x,Ql) = T]dAdx"-''xdy{x,Ql) , 

xASea{x,Ql) = risAsx""'xSea{x,Ql) , 

xAG{x,Ql) = r]gAgx''^{l-x)''^xG{x,Ql) (35) 

where on R.H.S. of ( pS]) we use the MRS(A') unpolarized densities. The normalization 
factors Af are determined in such a way as to ensure that the first moments of the 
polarized densities are given by r]f. Since the present polarized experiments do not 
allow for a flavour decomposition of the sea, SU(3) symmetry of the sea-quark densities 
is assumed 

Au{x, Ql) = Ad{x, Ql) = As{x, Ql) = Aq{x, Ql) . (36) 

Following this assumption ASea{x, Ql) = 6Aq{x, Qq) and rjs = Qrjq , where rjq is 
the first moment of Aq . Note also that the very small charm contribution to the 
unpolarized sea on R.H.S. of (^) is ignored in our analysis. 
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In this approach, the first moments of the valence quark densities rju and rjd are 
obtained directly from the sum rules (H) and (0) 

rju = 0.918, rid = -0.339 . (37) 

The first moment of the polarized sea ris is fixed from the measured value of Fi. 
More details are discussed in the next section. The rest of the parameters 

{ttu, ad, as ,r/g , % ,bg} , (38) 

have to be determined from the fit to the A^{x, Q"^) data. 

Finally it should be noted that using a set of polarized parton densities like (p5D, 
A(it,(a;, Qq) and ^q{x,Ql) will turn out to be negative in the whole x region in 
contrast to the BBS model where these quantities become positive at large x. 



5. Results of Analysis 

In this section we present the results of our fits to the world A^{x, Q^) data: EMC 
proton data 0, SLAC E142 neutron data 0], SLAG E143 proton and deuteron data 
[5 - 7], SMC proton data [Q and the SMC deuteron data |jl2| which are combined data 
from the 1992 |jlO|, 1994 [|ll|] and 1995 runs. The data used (203 experimental points) 
cover the following kinematic region: 

0.004 <x< 0.75, 1 < Q2 < 72 GeV^ . (39) 

We have chosen Qq = 4 GeV^. In all fits only statistical errors are taken into 
account. "Higher twist" corrections are not included in the present study. As already 
discussed above, in the approach used their effect is expected to be negligible. 



(a) BBS model 



A comparison of our results (solid curves) with the data on A^{x, Q"^) is shown in 
Figs. 2a - 2f. Our NLO results for g^^x^Q"^) are illustrated in Fig. 7. The minimum 
of the functional is achieved at = 232.7 and x^/DOF = 232.7/194 = 1.20. 
The values of the free input parton parameters ( p4D corresponding to this value 
are 

S„ = -3.010 ±0.156, C„ = 2.143 ±0.137, = -2.065 ± 0.148 , 

Cd = 1.689 ± 0.227, = 0.334 ± 0.044, = -0.292 ± 0.042 , 

Bg = -0.339 ± 0.454, a, = 1.313 ± 0.056, ag = 1.233 ± 0.073 . (40) 
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The rest of the parameters Au, Ad, Bd, Dd, Ag, Bg and Ag are determined 
by the constraints (0) and (0), the sum rules (|) and (H) and the momentum-energy 
sum rule, giving 

A^ = 3.088, Ad = 0.343, Bd = -0.265, Dd = -1.610 , 

As = 0.001, Bs = 0.041, Ag = 1.019 . (41) 

It is seen from Figs. 2a - 2f that the NLO QCD predictions using the BBS model 
for the input parton distributions are in a good agreement with the presently available 
data on A^{x,Q^) , as well as with the corresponding g^{x,Q^) data (see Fig. 7). 

For the first moments of the polarized flavour and fiavour-singlet distributions at 
Qo = 4 GeV"^ we obtain 

Au + Au = 0.839, Ad + Ad = -0.405, 2As = -0.079 , (42) 

and 

AS = ao = 0.342 . (43) 
These values yield for the quantity r^(Q^) , the first moment of Qi , 

r?(5 GeV'^) = 0.146, 1^(5 GeV"^) = -0.047, rf(5 GeV^) = 0.046 , (44) 

which are in a good agreement with their experimental values 

r?(5 GeV^) = 0.141 ± 0.011 {All proton data) , 

r^(5 GeV^) = 0.039 ± 0.004 {All deuteron data) . (45) 

The experimental values of and Tf in (|45|) have been determined in |^ and 
by combining SMC and SLAG E143 results on Ai and Af, respectively. 
Our result for the axial charge ao 

ao(5 GeV^) = AS(5 GeV^) = 0.341 (46) 

is also in a good agreement with its experimental value ao(5 GeV"^) = 0.29 ± 0.06 
recently determined from the combined analysis of all proton, neutron and deuteron 
data. 

We obtain for the small x behaviour of the input sea quark and gluon distributions 

xAq{x) ~ x°-^\ xAG{x) ~ . (47) 

This result confirms that of Gehrmann and Stirling who used in their NLO anal- 
ysis the same renormalization scheme and the same Ql , but different input parton 
densities. 
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It is well known that in the unpolarized case the gluon distribution can not be well 
determined by the fit to the data on nucleon structure functions alone. The situation in 



the polarized DIS is even worse. In particular, the fact that the parameter Bg in (^0|) 
is not well determined reflects this uncertainty in the extraction of AG{x, Q"^) from 
the data. We obtain for the first moment of AG(x,Q2) : AGi(4 GeV"^) = 0.447 , 
which does not coincide with the mean values of this quantity given in most of the 
theoretical analyses. This result for AGi{Ql) not surprising. It should be noted 
that in the BBS model one can show that AGi{Ql) constrained by 

AGM) < ^G2{Ql) - G^iQl) , (48) 

where G2{Ql) and G^IQq) are the second and third moments of the unpolarized 
densities. Bearing in mind that from the unpolarized data G2(Qo) ~ 0-^5 one has 
that AGi{Qq) should be smaller than 0.9. This fact could be critical for the validity 
of the BBS model if more precise data are available in the future. 



The BBS polarized x{Af{x) + Af{x)) distributions at Ql = 4 GeV^ are shown 
in Fig. 3. In Fig. 4 we compare our results for the BBS unpolarized x{f{x) + f{x)) 
input densities (solid curves) with the MRS(A') set of unpolarized parton distribu- 
tions (dashed curves) at the same Qo = 4 GeV^. The difference between the BBS 
unpolarized input parton densities and the MRS (A') ones is somewhat greater than 
the usual difference between the various sets of unpolarized parton distributions used 
in the literature. The difference for x{u + u) and x{d + d) quarks increases up to 20% 
depending on x, while e.g., for MRS(DO') and MRS(A') set of input valence quark 
densities it is typically 6-8% in the x range: 0.01 < x < 0.50 and amounts to 15-20% 
for the range: 0.5 < x < 0.7. As one can see from Fig. 4 the agreement for the sea 
quarks and the gluons is better, of the same quality as for the well known different sets 
of unpolarized input partons. 



(b) Models based on the unpolarized parton densities 



Let us continue now with discussion of our results of the fit to the data using for 
the input polarized parton densities the set (p5|). Such a set of input partons had been 
used in ||2^, but starting at very small Ql = 0.34 GeV'^. 

Unlike the BBS model the data do not allow to determine r]s, the first moment of 
the polarized sea, in proper way if the parametrization (p5D is used. In order words, 
if ris is taken to be a free parameter, its value determined from the fit to A^{x, Q^) 
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data, does not agree with the experimental value of r^s/S = —0.10 ± 0.02 [Q. That 
is why we fix r]s from the measured value of rf(5 GeV"^) = 0.141 ± 0.011 to be 



Vs 



-0.290 at Ql = 4:GeV^. 



As already mentioned above it is impossible to determine accurately the form of the 
polarized gluon density from these data alone and therefore additional constraints have 
to be applied to the gluonic input parameters. As in |^ we have used the assumption 



ttg = Qs which defines the behaviour of AG'(a;, Q^) at small x. It turns out that even 
in that case it is not possible to determine by the fit the value of the parameter hg in 
( p5| ) which controls the behaviour of AG(a;, Q^) at large x. For that reason the fit to 
the data was performed at different fixed values of hg in the range < 6^ < 7. Further, 
if ad is left as a free parameter, its best-fit value is negative independently of the value 
of hg in the above range. This fact, presumably refiects the circumstance discussed 
in Section 2, namely that from the data in the small x region gi — gi > gl + gi , 
which probably induces such a behaviour in Adv{x,Ql) . For this negative value of 
ad'-ad^ -0.13 

= 0.33M,.- 

is greater than 1 at very small x : x < 0.3.10"^ and the positivity condition (P) is 
broken. In addition, the corresponding value of determined by this fit is such that 
the positivity condition (||) for ASea{x) is violated too, but at large x : x > 0.80. 
That is why we have taken = , the smallest value of which guarantees positivity 
for at all x. The change of from x^i^d < 0) = 215 to x^i^d = 0) = 220 is 
negligible. 

The results of the fit to Ai{x, Q"^) data are presented in Table 1. 

Table 1. The results of the NLO QCD fit to the world data using the set 
( p5D for the input polarized partons ( = 0). 





hg = 


hg=5 


hg = 7 


xVdof 


221.2/200 


220.5/200 


219.5/200 




0.196 ± 0.018 


0.178 ± 0.027 


0.169 ± 0.025 


as 


0.892 ± 0.085 


0.706 ± 0.155 


0.694 ± 0.112 


Vg 


-0.14 ± 0.26 


1.07 ± 0.70 


1.29 ± 0.53 



Although only three free parameters {a^, a^, rjg} have been used a very good 
description of the data is achieved. This is illustrated in Figs. 2a - 2f (dashed curves). 
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which show the NLO description of the various measurements. In Fig. 7 we 

compare our NLO results for g^i^x^Q"^) with the SMC data. The input polarized 
parton distributions ( ^5|) are shown in Fig. 5. The value of x^/DOF is 220.5/200 
for bg = 5 and practically does not depend on bg (see the Table). This value is 
better than that one in the case of the BBS model. The mean value of rig , the first 
moment of the polarized gluon density, is sensitive to the value of bg : rig increases 
with increasing bg. The values of rig at bg > 5 are greater than 1 and are in 
agreement with those determined by the NLO analysis in The sets of polarized 
quark densities corresponding to gluons with bg = 5 and bg = 7 are practically the 
same. 

We have found also that values for rig smaller than 1, and even a small negative 
value rig = —0.14 ± 0.26 corresponding to bg = , are not excluded by the data. In 
the last case, however. 



\^Sea{x)\ _ ^ n.892 



Sea{x) 

and positivity is violated in the large x region: x > 0.88. 



In Fig. 6 a comparison between the BBS input polarized distributions ( ^0|) and the 
polarized densities (|35|) is shown. It is seen that while the distributions x{Au + An) 
and 2xAs = 2xAq are very similar for these two parametrizations, the polarized dis- 



tributions x{Ad + Ad) and xAG corresponding to ( pOf ) and ( p5| ) differ considerably, 
x{Ad + Ad) for x > 0.35 and xAG for all x. The polarized gluon distribution enters 
gi{x, Q^) at next-to-leading order and bearing in mind the accuracy of the present data 
these different gluon contributions can not be distinguished. However, the difference 
between x{Ad + Ad){x, Ql) leads to a considerably different behaviour of Af{x, Q"^) 
in the kinematic region: x > 0.35, ~ 5 — 10 GeV^, and allows a better fit to the 
SLAG E143 data on Af{x, Q^) in this region in the case of the parametrization ( |35D 
for the input polarized parton densities. The difference for x{Ad + Ad) at large x is 
a consequence of the fact that the BBS distributions are forced to satisfy (pSf) as s ^ 1. 



5. Conclusion 

We have performed a next-to leading order QCD analysis {MS scheme) of the 
world data on polarized deep inelastic lepton-nucleon scattering. The QCD predictions 
have been confronted with the directly measured virtual photon-nucleon asymmetry 
A^ [x, Q"^) rather than with the polarized structure function gi{x, Q"^) derived from 
the data by different additional procedures. Different parametrizations for the input 
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polarized parton densities have been examined: the Brodsky, Burkardt and Schmidt 
model and a model based on unpolarized parton densities. In the BBS model the 
positivity constraints for the unpolarized parton densities are automatically valid. In 
the second one it is easy to control them. In both cases a good fit to the data is 
achieved. The description of the data is slightly better if for the input polarized parton 
distributions the second parametrization is used. 

A distinctive feature of the BBS model is that the helicity-dependent input parton 
distributions f±{x,Ql) (rather than Af{x,Ql) and f{x,Ql) ) are parameterized 
so that the polarized, as well as the unpolarized data have to be fitted by the same 
parameters (|3^). This feature creates a serious challenge to the model. We have shown 
that the BBS model is compatible with the present DIS data although the agreement 
with the MRS (A') parametrization of the unpolarized parton distributions at = 
4 GeV^ is somewhat worse than the usual level of agreement between the different 
sets of the input unpolarized parton densities usually used in the literature. Our result 
for the axial charge ao(5 GeV'^) = AS(5 GeV'^) = 0.341 , obtained in the BBS model 
is in a good agreement with its experimental value ao(5 GeV"^) = 0.29 ±0.06 recently 
determined from a combined analysis of all proton, neutron and deuteron data, while 
using the second parametrization for the input polarized densities this quantity is not 
well determined from the fit to the data. This results from the fact that the sea quark 
distributions are still largely undetermined if the second kind of parametrization is 
used. 

The present polarized data do not allow a precise determination of the shape of the 
polarized gluon density. It follows from our fits that values of AG'i(4 GeV"^) , the first 
moment of the polarized gluon density, both greater and smaller than 1, are possible. 
Negative values of AGi are not excluded either. 

Despite the great progress of the past few years it is clear that in order to test 
precisely the spin properties of QCD more accurate DIS polarized data at fixed x and 
various are needed. In addition, charged current data will be very important for 
a precise determination of the polarized parton densities and especially, for a precise 
fiavour decomposition of the polarized quark sea. Finally, a direct measurement of 
AG'(a;,Q^) in processes such as J/ip production in lepton-hadron scattering with a 
polarized beam will answer the important question about the magnitude and the sign 
of AGi. 
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Figure Captions 



Fig. 1. Comparison of the SMC data ||^, 
value = ig GeV^ . 

Fig. 2a-2f. Comparison of our NLO results for A^{x,Q'^) with the present data. The 
solid and dashed curves are the best fits to the data corresponding to the BBS model 
(|30| ) and parametrization (|35D of input polarized parton distributions, respectively. 

Fig. 3. Next-to-leading order input polarized parton distributions at = 4 GeV^ 
in the BBS model (xAm = xAu + x^u and xArf = xArf + xAd). 

Fig. 4a-b. Comparison between BBS unpolarized parton distributions (solid curves) 
and MRS(A ') set of unpolarized parton densities (dashed curves) at = 4 GeV^ {xu = 
xu + xu and xd = xdu + xd) . 

Fig. 5. Next-to-leading order input polarized parton distributions at = 4 GeV^ 
corresponding to parametrization ( PBj ) in the text (a^ = 0, bg = 7). 

Fig. 6. Comparison between the input polarized parton distributions in the BBS 
model (solid curves) and parametrization ( |35D (dashed curves), xu = xu + xu, xd = 
xdu + xd and 2xAs = 2xAq. 

Fig. 7. Comparison of our NLO results for g^{x,Q^) with SMC data [g [g at 
the measured Q^. The solid curves correspond to BBS model and dashed ones to 
parametrization (^) of the input polarized parton densities. 



13] on (^i ± gi as function of x at mean 
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